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Abstract: We study a static black hole localized on a brane in the Randall-Sundrum (RS) 
II braneworld scenario. To solve this problem numerically, we develop a code having the 
almost 4th-order accuracy. This code derives the highly accurate result for the case where 
the brane tension is zero, i.e., the spherically symmetric case. However, a nonsystematic 
error is detected in the cases where the brane tension is nonzero. This error is irremovable 
by any systematic methods such as increasing the resolution, setting the outer boundary at 
more distant location, or improving the convergence of the numerical relaxation. We discuss 
the possible origins for the nonsystematic error, and conclude that our result is naturally 
interpreted as the evidence for the nonexistence of solutions to this setup, although an 
"approximate" solution exists for sufficiently small brane tension. We discuss the possibility 
that the black holes produced on a brane may be unstable and lead to two interesting 
consequences: the event horizon pinch and the brane pinch. 
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1. Introduction 

The Randall-Sundrum (RS) braneworld scenarios Q |2) have been attracting a lot of at- 
tentions since their appearance in 1999. In these scenarios, our world is a 5-dimensional 
anti-de-Sitter (AdS) spacetime with a negative cosmological constant A. In the RS I sce- 
nario [fflj, there are two branes with positive and negative tension, and our 3-dimensional 
space is the negative tension brane. In this scenario, the Planck energy could be O(TeV) 
and the black hole production might happen at the Large Hadron Collider (LHC) |], |5|. 
In the RS II scenario there is only one brane with positive tension and the extra- 
dimension is infinitely extended. By the bulk curvature effect, the 4-dimensional gravity 
is successfully realized on the brane || [jj. Although the hierarchy problem is not solved, 
this scenario is also interesting since it gives the new mechanism of the compactification of 
the extra dimension. 

In this paper, we study a static black hole localized on the brane in the RS II scenario. 
There are at least three motivations for this study. The first motivation is related to 
the astrophysics in this scenario. The final state of the gravitational collapse of a star is 
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expected to be different from the 4-dimensional black hole, since the event horizon extends 
into the bulk. By studying black hole solutions, we might obtain observable features of 
braneworld black holes or constraints for the RS models. The second motivation is the 
black hole production at the LHC Although the study of the black holes in 

the RS II scenario is not directly related to this phenomena, the obtained results can be 
translated into the RS I case by changing the sign of the bulk curvature scale I if the effect 
of the existence of the positive tension brane is neglected. Therefore, we would obtain 
some indication for the effect of the bulk curvature on the black hole physics that might 
be observed at the LHC. 

The third motivation is the AdS/CFT correspondence, which conjectures that the 
gravitational theory in the AdS spacetime is dual to the conformal field theory (CFT) 
on the boundary of that spacetime. It was suggested that the AdS/CFT correspondence 
hold also for the RS II scenario [||, ||, 10]. Based on the AdS/CFT correspondence, it was 
conjectured that no solution of a large black hole on a brane exists pi] , |L2]| . Their argument 
is that a 4-dimensional black hole with quantum fields is dual to a classical 5-dimensional 
black hole on a brane, and the dual phenomena of evaporation of a 4-dimensional black hole 
is expected to be escape of a 5-dimensional black hole into the bulk. However, there are 
different opinions on how to apply the AdS/CFT correspondence and thus on the expected 
black hole solution space [||, U, |l|, |l|, 0, HI. The explicit study of the black hole 
solution would shed a new light on this issue. 

There were many attempts to obtain solutions of a black hole localized on a brane. 
The exact solution of a static black hole intersecting a 2-brane in a 4-dimensional space- 
time was found in [19, 20 1 using the so-called C-metric |21]. The brane is given as some 
(2+l)-dimensional slice in the spacetime of the C-metric, and the geometry of the brane is 
asymptotically a cone. In contrast to the 4-dimensional case, finding solutions of a black 
hole on a brane turned out to be very difficult in the 5-dimensional case. For several years, 
no successful discovery of the solution was not reported, although several efforts were made 



m, 23, E3, 25, 26, 27, 



After that, Kudoh et al. [29] reported a numerical study of a static black hole on 
a brane. In their formulation, the problem was reduced to elliptic equations for metric 
functions with appropriate boundary conditions. They solved this problem by a relaxation 
method. For the cases where the horizon radius is much smaller than the bulk curvature 
scale (L := l/ph 3> 1), they realized the sufficient convergence. As the mass of the black 
hole becomes large, the convergence became worse and the error grew. For this reason, 
they only showed the results for 3 < L < 500. Since their study, no calculation has been 
performed for the case L < 1. The result of 2£] might be interpreted as an evidence for 
the existence of solutions of a black hole on a brane with a small mass. However, the 
growth of the error with (1/L) also can be regarded as the evidence for nonexistence of 
such solutions. Therefore, we have to be careful in interpreting their results, and this study 
has to be reexamined. 

The solution of a static black hole on a brane was studied also by the perturbative 
method by Karasik et al. [30, 31] . They studied the situation where the black hole radius 
is much smaller than the bulk curvature scale using the matching method, and derived 
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the solution that is at least C 1 on the horizon. In Appendix A, we summarize the part of 
their results that is closely related to the discussions in this paper. As pointed out there 



(and also in the original paper [31 1), there remains a possibility that the horizon becomes 
singular at the higher-order perturbations. Therefore, this study also cannot be regarded 
as the rigorous evidence for the existence of the solution. Another related study is the 



perturbative study of a higher-dimensional version of the C-metric by Kodama [32]. Using 
the gauge invariant formulation for the perturbation of a higher-dimensional Schwarzschild 
black hole |33|, [34]], he succeeded in deriving the solution of a 5-dimensional black hole 
accelerated by a string with uniform tension. However, in contrast to the 4-dimensional 
case, there was no (3+l)-dimensional slice that realizes the junction condition on the brane 
in that spacetime. 

As summarized above, there is no rigorous answer to the solution space of a static black 
hole on a brane. We have to examine the (non) existence of the solution, and if it exists, 
the range of L < 1 has to be explored. Although we basically follow the formulation by 



Kudoh et al. [29], a significantly improved numerical method is developed. Specifically, we 
focus our attention to the error analysis taking account of the possibility that no solution 
exists. It is usually said that numerical calculations cannot prove nonexistence of solutions 
rigorously. However, since our code has the sufficient accuracy, it is possible to obtain a 
strong indication for (non)existence of the solution by such an error analysis. In this way, 
we obtain the result that at least raises strong doubts for the existence of solutions of a 
black hole on a brane. 

This paper is organized as follows. In Sec. ^ we explain the setup of the problem, which 
is a short review of the corresponding part of Ref. P9fl . In Sec. our numerical method 
is explained taking attention to the difference from the code of Ref. [p9| . In Sec. ^j, we 
show the results for the case where the brane tension is zero, i.e. the spherically symmetric 
case. This proves the correctness and high accuracy of our code. After discussing the 
features that the numerical solutions have to satisfy, we introduce "systematic error" and 
"nonsystematic error", which are the key notions in order to interpret our results for the 
cases where the brane tension is nonzero. In Sec. ^J, the numerical results for the cases of a 
nonzero tension brane are shown taking attention to the numerical errors, and the existence 
of the nonsystematic error is proved. We also point out that our results are consistent with 
the perturbative study if the numerical errors are ignored. Then, we discuss the 

interpretation of our results. Our results are naturally interpreted as the evidence for the 
nonexistence of solutions of a black hole on a brane. Section || is devoted to summary and 
discussion. Assuming the nonexistence of solutions of a static black hole on a brane, we 
discuss possible phenomena that could happen after the black hole formation on a brane. 
In Appendix |A|, we briefly summarize the results of the perturbative study |3C, 31 1 of a 
small black hole on a brane, taking attention to the part that is closely related to this 
paper. 

2. Setup 

In this section, we explain the setup to study a static black hole localized on a brane, which 
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is spherically symmetric on the brane (in the 4-dimensional sense) and axisymmetric in the 
bulk spacetime (in the 5-dimensional sense). Since we basically follow the formulation of 



Kudoh et at, this is a brief review of the corresponding part of Ref. [29]. 
We assume the following metric ansatz: 

& 2 = ^ \_ T 2 dt 2 + e 2R {dr 2 + + ^20^2] ^ (2J) 

where dQ 2 := d6 2 + sin 2 9d(p 2 . Here, £ is related to the bulk cosmological constant as A = 
—6/£ 2 . The functions T, R, and C depend only on z and r. Each hypersurface z = const, is 
a 4-dimensional spherically symmetric spacetime where r is the radial coordinate, r = is 
the symmetry axis in the 5-dimensional point of view. If we set T = 1 and R = C = 0, the 
spacetime is the AdS spacetime in the Poincare coordinates. In the metric fl2.1| ), the fact 
that the timelike killing vector is hypersurface orthogonal is used following the standard 
definition of a static spacetime (e.g. Ref. (3j|). The fact that a 2-dimensional surface is 
conformally flat is also used in (|2.1| ). For these reasons, the assumed form of the metric 
( |2.l| ) is sufficiently general for our setup. 

We introduce the coordinates p and x by 

r = psmx'i z = £ + pcosx- (2-2) 
Then, T, R, and C are the functions of p and x- The locations of the event horizon and the 



brane are assumed to be p = ph and x = 7r /2, respectively. Note that it was shown in [29 
that this assumption is possible in general using the remaining gauge degree of freedom in 
the metric. Although there are the two parameters £ and ph in this system, the spacetime 
is basically specified by the dimensionless parameter 

L := £/p h , (2.3) 

since fixing L and changing £ or p^ corresponds to just changing the unit of the length. In 
the numerical calculation, we put the event horizon at ph = 1- In this choice of the length 
unit, we have L = £. In this paper, we basically use (1/-L) as the parameter to specify the 
spacetime rather than L. 

2.1 Equations 

The equations to be solved are the Einstein equations of a vacuum spacetime with a negative 
cosmological constant, 

Qy.v ■= R^u ~ (2/3)A^ = 0. (2.4) 

From Q\ = 0, Q\ - G p p - Q% + 2Q e e = and Q% = 0, the elliptic equations for T, R, and C 
are obtained: 

7 2rr , n („ , 2^ \\„ , 2 , „ , 2p 



V Z T + 2 ( C p + — - - ) T p + ( cot X + C, X + -f sin X ) T t> 



2 4 / K( 2 \ 

+ — (sin xC, x - cos X pC, p ) T + f 1 + — e 2R J T = 0; (2.5) 
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+ C p (2C p + I - + ^ ( 2 C x + 4 cot x + ^ sinx) = 0. (2.7) 

Here, V 2 := d 2 + d p /p + <9 2 /p 2 . There are two other equations derived by Q p q = and 
Q\ -Qp + Qx+ ig 6 e = 0, respectively: 
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0. (2.9) 



It was shown in that a solution of Eqs. fl2.5|) -(2.7) automatically satisfies Eqs. (J 
and ( |2.9D provided that Eq. is satisfied on the boundary and Eq. ( |2.9|) is satisfied at 
one point. Therefore, the latter two equations ( |2.8| ) and ( |2 . 9| ) are the constraint equations, 
and the main equations are the first three equations ( j2.5| )-( ^77| ). 
In the limit (1/L) := Ph/£ — ► 0, the solution is 



P 2 ~l 

P 2 + r 



R = Cn = log 1 + 



I 



(2.10) 



which gives the metric of the Schwarzschild black hole in the isotropic coordinates. This fact 
raised the expectation that the Schwarzschild- like black hole solution exists for (1/L) <C 1 



in the previous studies [£9], |3fj, |31 1 . 
2.2 Boundary conditions 

There are four boundaries of the computation domain: the outer boundary p = p ou t, the 
brane X = 7r /2, the symmetry axis x = 0, and the horizon p = ph. Here we summarize the 
boundary conditions for T, R, and C on the four boundaries one by one. 
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At the distant region, the spacetime should be reduced to the AdS spacetime. There- 
fore, we impose 

T = l; C = 0; R = 0, (2.11) 
at p = pout for a sufficiently large p out . Strictly speaking, we have to impose the conditions 



( 2.11 ) at p = oo. Since we adopt the finite value of p ou t for numerical convenience, we have 
to keep in mind the existence of an error from the finiteness of p ou t (say, the "boundary 
error") and check the dependence of the numerical errors on the value of p ut- 

On the brane x = 7r /2, the boundary conditions are given by Israel's junction condition 

K^u = (2.12) 

where K^ v and 7^ are the extrinsic curvature and the induced metric on the brane, 
respectively. This condition implies that the brane has the constant tension a = 3/ (^ttG^I). 
The condition fl2.12|) is expressed as 

^ = R x = C x = P £ (e R - 1) (2.13) 

in terms of the metric functions, and the constraint equation ( |2.8| ) is manifestly satisfied 
on the brane by these boundary conditions. 

At the symmetry axis x = 0, we impose the regularity conditions 

T x = R x = C x = 0, (2.14) 

and 

R = C. (2.15) 



The second condition ( 2.15 ) is implied by the regularity of the second terms in Eqs. ( 



and (2.7). Under the conditions ( 2,15| ) and ( |2,14 ), the first constraint equation ( |2.8[ ) is 



trivially satisfied at \ = 0. Substituting Eqs. (|1|) and (|2~li ) into Eqs. (^5|) -(^) 
we obtain the formulas for T xx , R, xx , and C tXX . The second constraint equation 



evaluated at x = is consistent with these formulas. In our numerical calculation, all 



these conditions are imposed as explained in Sec. 3.1 



On the horizon p = ph, the timelike killing vector should become null, and thus 

T = 0. (2.16) 
Then, the regularity of Eqs (|2.6|)-( p^9[ ) implies 

C P = --; (2.17) 

zp 

T ,PX = R ,X T ,P> ( 2 - 18 ) 



and 



3/ 

2R tP + C„ = . (2.19) 

zp 



- 6 - 





Our code 


The code of Ref. g|] 


radial coordinate 


x := log p 


p (geometric progression) 


angular coordinate 


X 


C ■= x 2 


scheme 


almost 4th-order 


2nd-order 


grid size 


Ax = 0.025-0.05 


A Pl = 0.049 


grid error 


~ 10 _7 -10~ 5 , 


~ 10~ 2 


outer boundary 


4 < Xout < 7 


Pout = 85 


variables 


T,X,Y 


T, R,C 



Table 1: Comparison between our code and the code by Kudoh et al. [ p9[ . 



Here, the condition ( |2.18| ) is equivalent to the zeroth law of the black hole thermodynamics, 
i.e. the constancy of the surface gravity 



-Rr 



(2.20) 



on the horizon, and the condition ( [2.19 ) is equivalent to the zero expansion of the null 
geodesic congruence on the horizon. Here, it has to be pointed out that the four condi- 
tions exist for the three variables. If a regular solution exists, it has to satisfy all these 
four boundary conditions on the horizon. But in the numerical computation, only three 
conditions can completely fix the solution to Eqs. fl2,5| )-( |2~7r] ). The three conditions ( 2,16j ), 
fl2~T7D , and (|2~T9| ) are used in our code, since the use of Q2.18D leads to a numerical in- 
stability. Therefore, the constancy of the surface gravity k is not explicitly imposed in 
our calculation, and whether it is satisfied (within the expected numerical error) has to be 
checked after generating a solution. 



3. Numerical method 

In this section, we explain how to solve the problem numerically. Since our method is 
significantly improved compared to the previous study by Kudoh et al. p9[ , we focus 
attention to the differences between the two codes, as summarized in Table [l|. 

3.1 Coordinates, variables, and scheme 

In the code of Ref. |29| ], the nonuniform angular coordinate £ = X 2 was adopted, because 
the second terms proportional to 1/ sin 2 ^ in Eqs. ( |2.6[ ) and (2.7) cause a severe numerical 



instability and the introduction of £ makes the problem more tractable. But in this coor- 
dinate, the grid size measured in x ls n °t so small in the neighborhood of the symmetry 



axis. For the grid number 100 in |29], the grid next to the symmetry axis is located at 
X = vr/20 ~ 0.16, and this leads to a large numerical error. For this reason, we adopt the 
uniform angular coordinate x m our code. In this case, some treatment is required in order 
to realize the numerical stability, as explained in the next subsection. 

The choice of the radial coordinate is also different. The authors of [^] adopted the 
method of the geometric progression, in which the grid size becomes larger as the value of 
p is increased. They put the outer boundary at p out = 85 and used the grid number 1000, 
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where the ratio of the last grid size to the first grid size was set as Ap max /Api ~ 2.71 
with Api = 0.049. For this choice, the expected grid error is ~ (Api) 2 ~ 0.24%, and it 
is enhanced to ~ 1% in the neighborhood of the symmetry axis. This is not so good as a 
numerical calculation. In our code, we use the nonuniform radial coordinate x defined by 

x:=\og{p/p h ). (3.1) 

This means that the grid size measured in p is proportional to p. However, this coordinate 
choice does not cause a significant numerical error, because the error becomes large only 
if the functions change rapidly, while the functions T, R, and C asymptote to the constant 
values at the distant region. To be more precise, the functions that decay as ~ 1/p 
behaves as ~ exp(— x) in the coordinate x, and such an exponential decay is numerically 
tractable. For this reason, we obtain the sufficient numerical accuracy with relatively 
small grid numbers. In our calculation, the outer boundary is located at x on t = 10 (i.e. 
/5 out — 22026.5) for the case of a zero tension brane (Sec. Q). In the cases of a nonzero 
tension brane, W6 chcinge the value of #out 

around 5 (i.e. p OVLt — 148.4) and observe the 
dependence of the error on x out (Sees. O] and |5.2|). 



Let us turn to the choice of the variables. In Ref. [29], the variables T, R, and C were 



solved in the numerical calculation. In our code, we introduce 

X:=R + C; Y:=R-C, (3.2) 

and solve T, X, and Y. There are two reasons for this choice. The first reason is related 
to the numerical stability. From Eqs. ( [2.6D and ( |2.7j ), we derive the equations of the form 
V 2 X = ••• and V 2 Y = In these equations, the term proportional to l/sin 2 x is 

included only in the latter equation, and thus it becomes easier to handle the instability 
caused by this term. The second reason is that the nonsystematic error in the case of a 
nonzero tension brane appears as an unnatural jump in the value of Y in the neighborhood 
of the axis. 

To summarize, the variables T, X, and Y are solved using the coordinates x and 
X in our code. We adopt the finite difference method, where the grids are located at 
x = IAx (I = 0, ...,/max) and x = J^X {J = 0, J max ). Here, Ax := x out /I max and 
Ax '■= (7r/2)/J m ax> and we keep the ratio Ax/ Ax = ir/2. We adopt the 4th-order accuracy 
scheme except at the grids in the neighborhood of the three boundaries x = 7r /2, x = 0, 
and x = x ou t; where the 3rd-order accuracy scheme is adopted. Therefore, our numerical 
code has the almost 4th-order accuracy, and therefore we can obtain more accurate results 
with smaller grid numbers compared to the case that the 2nd-order accuracy scheme is 
adopted. As proved in next section, the error by the finiteness of the grid size (say, the 
"grid error") is < 10 -5 for Ax = 0.05 and ~ 10~ 7 for Ax = 0.025. Compared to the grid 
error ~ 1% in Ref. (2^], the numerical accuracy is greatly improved. 

We briefly mention how to impose the boundary conditions at the symmetry axis. 
There, we have the regularity conditions Y = 0, T x = X x = Y tX = 0, and formulas for 
T xx , X )XX , and Y tXX . It is possible to write down the finite difference equations for T, X, 
and Y at J = (i.e. x = 0) that include the conditions for the first and second order 
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derivatives. We adopt those finite difference equations for T and X in order to determine 
the values of T and X at J = 0. As for Y, we impose Y = at J = 0, and use the remaining 
finite difference equation for Y to determine the value of Y at J = 1 (i.e. x = A\0- This 
is better than using the finite difference equation at J = 1 that has the term proportional 
to 1/ sin 2 x- 

3.2 The numerical relaxation 

In order to solve the finite difference equations, we adopt the relaxation method, in which 
one prepares an initial surface and makes it converge to the solution iteratively. Our code is 
based on the successive-over-relaxation (SOR) method. Namely, we calculate the difference 
from the finite difference equations AT(jj), Al^^j, and AY^j), and determine the next 
surface using the formulas 

r a7) )=r (^)+^AT (/ , J); 

X M = X W + ™xAX M , (3.3) 
Yffi =Y {I ,j ) +w Y AY {I>J) . 

In the usual SOR method, wt, wx, and wy are called the acceleration parameters and 
are chosen to be a value between 1 and 2. In our case, in order to realize the stability, we 
choose 

wt = wx = 0-8; wy = 0.08 x cos Q [(x - x out )Tr/2] sin x, (3.4) 

where a is chosen appropriately depending on the situation. The reason for the complicated 
functional form of wy is as follows. The numerical instability tends to happen at the distant 
region or near the symmetry axis. In those regions, the value of wy is very small and this 
makes the convergence of the surface very slow. In this way, we can avoid the instability 
caused by the term proportional to 1/ sin 2 x- 

In order to evaluate the degree of the convergence, we introduce the parameters 

f EjA^ol , _ e|ax (j , j) | j _ E \ ay (i,j)\ 

€T = Ft P x ~ v- I y P ey ~ P^ F • ^• 5 ' ) 

Z^| J (/,J)| Z^| A (7,J)| JmaxJmax 

In the formula of ey, the normalization factor was not adopted as Yl Y (i,J) because the 
analytic solution is Y = in the case of a zero tension brane. Typically, the value of ey is 
the largest among the three convergence parameters. We truncated the relaxation process 
when 

max(eT, ex, ey) < eo (3.6) 

is achieved. In the cases of a nonzero tension brane, the value of eo is taken as eo = 10~ 10 . 
The error from this truncation (say, the "relaxation error") is ~ 10~ 6 . In the case of a 
zero tension brane in the next section, we adopt eo = 10~ 12 where the relaxation error is 
~ 10~ 8 , because we have to make the relaxation error sufficiently smaller than the grid 
error in order to prove the appropriate convergence of the numerical solutions. 
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Figure 1: The relation between the grid size Ax and the error 5 in the case of a zero tension brane 
in the log scales. The numerical data shows the almost 4th-order convergence. 

4. Results for a zero tension brane 

In this section, we show the results for the case of a zero tension brane, i.e. the spherically 
symmetric case. In this case, (1/L) = and the analytic solution ( 2.1CQ (i.e., Tq, Rq, and 



Co) exists. A comparison of the numerical data with the analytic solution gives a useful 
check for the correctness of our code. Furthermore, we can learn the features that the 
numerical solutions have to satisfy when the solutions are actually present. 

In these calculations, we set p/£ = in our code. Except this, we do not explicitly 
impose the spherical symmetry. As a result, the numerical solution slightly depends on x 
because of the numerical errors. It is worth pointing out that the sufficient convergence 
was realized with the initial surface that depends on x, in spite of the presence of the term 
proportional to 1/ sin 2 %■ This proves the sufficient stability of our code. We calculated T, 
X and Y for various locations of the outer boundary x out and grid sizes Ax. We changed 
the value of x ou t from 3 to 10 and obtained the natural result that the solution converges 
as Xout is increased. 

It is important to confirm the appropriate convergence of the numerical data. For this 
purpose, we change the grid size Ax from 0.025 to 0.2 for a fixed x ou t = 10. For this choice 
of x ou t, the boundary error is smaller than the grid error. Figure [l| shows the relation 
between Ax and the error defined by 

= E/,j(l^-?ol + l*-^ol + l^-^ol) 

Ei,j(\To\ + \Xo\ + \Y \) ' [ -> 

From this figure, the almost 4th-order convergence is confirmed. In fact, the slope of the 
curve in Fig. [l] is ~ 4.8, a bit larger than 4. This is interpreted as follows. Since we 
partly use the 3rd-order accuracy scheme, the error is expected to be between 0(Ax 3 ) and 
0(Ax 4 ). As Ax is decreased, the ratio of the number of the 3rd-order grids to the number 
of the 4th-order grids becomes smaller and this change in the ratio makes the error from 
the 3rd-order grids less effective. As a result, the slope of the curve becomes larger than 4 
in Fig. 
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Figure 2: The angular dependence of the error Sy(x) m Y f° r Ax = 0.1 (o), 0.05 (x), and 
0.025 (□). The error shows the almost 4th-order convergence and is almost invisible for the case 
Air = 0.025. The value of Sy(x) tends to decrease as the value of x is decreased. 

Let us discuss the angular dependence of the numerical error in Y, since it will be 
important in the case of a nonzero tension brane. Since Rq = Co and thus Yq = 0, the 
numerical value of Y itself represents the error. For this reason, we define 

5y(x) ■= max[y] x - min[y] x (4.2) 

as the error in Y characteristic to the angular direction x- Here, max[Y] x and min[Y] x 
mean the maximum and minimum values of Y for a given coordinate value x- Figure || 
shows the behavior of Sy(x) f°r the grid size Ax = 0.1, 0.05, and 0.025. The almost 4th- 
order convergence is again confirmed. Moreover, it is seen that 5y(x) is suppressed in the 
neighborhood of the symmetry axis. It is important to point out that this tendency is in 
contrast to the case of the 2nd-order code. In that case, the value of Y tends to increase 
as the value of x is decreased, and suddenly jumps to Y ~ at % = Ax because of the 
imposed boundary condition. Such an unnatural behavior does not appear in the 4th-order 
code, and it is one of the merits in using the 4th-order accuracy scheme. 

To summarize, we have proved that our code has the sufficient numerical stability and 
successfully reproduces the solution ( 2. 10] ) in the zero tension case. Our result shows that: 



(i) The numerical solution converges as the value of x ou t is increased; (ii) The numerical 
solution shows the almost 4th-order convergence; and (iii) The numerical error in Y is 
suppressed in the neighborhood of the symmetry axis x = 0. These three features are 
expected to be held also for the cases where the brane tension is nonzero as long as (1/L) 
is sufficiently small and regular solutions exist. 

Here, we introduce "systematic error" and "nonsystematic error, " which play impor- 
tant roles in the next section. The systematic error means the error that originates from 
the adopted approximation in the numerical method. The systematic errors were already 
introduced in this paper: the "boundary error" from the finiteness of the location of the 
outer boundary x ou t, the "grid error" from the finiteness of the grid sizes, and the "re- 
laxation error" from the truncation of the numerical relaxation. In our method, there is 
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Figure 3: The numerical solution of Y for (1/L) — 0.06 as a function of x and \- Here, x out = 5 
and Ax = 0.05 are adopted. The value of Y does not decay properly at the distant region. Besides, 
there is an unnatural jump in the value of Y between the grids J = 1 and 2 (i.e. x = A% an d 2A%). 



no other source of the numerical error. We can avoid the relaxation error by making the 
criterion of the truncation strict. The boundary error has to decrease as x out is increased, 
and the grid error has to decrease consistently with the adopted scheme as Ax is decreased. 
In other words, when the numerical convergence is sufficient, the obtained solution has to 
satisfy the features (i) and (ii) above. If we detect an error that does not satisfy the features 
(i) and (ii), we call it the nonsystematic error. If such a nonsystematic error is detected, 
we have to suspect that something is wrong. 



5. Results for a nonzero tension brane 

In this section, we show the numerical results for the case where the brane tension is 
nonzero. In Sees. |5.l| and |5.2|, the errors in Y and the surface gravity k are discussed in 



detail, respectively, and the existence of the nonsystematic error is shown. In Sec. 5.3 



we compare our results with the perturbative study [31] ignoring the numerical errors. In 



Sec. 5.4, the interpretation of our numerical results is discussed. 



5.1 The error in the variable Y 

There is a limitation in the values of (1/L) and x ou t for which numerical solutions can be 
obtained. For a fixed (1/L), the numerical relaxation is stable only for x ou t < x out (V-^) 
for some critical value x^. (1/L). The value of x^^\l/L) becomes smaller as (1/L) is 
increased and depends also on the resolution. When the numerical instability does not 
occur, we always realized the criterion ( |3.6| ) for the truncation of the relaxation. 

When both (1/L) and x ou t are not so large (e.g. (1/L) = 0.02 and x out = 5.0), the 
numerical solutions for T, X, and Y (apparently) look natural, and their gross features 
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Figure 4: Dependence of the error Sy in Y on the location of the outer boundary x ou t- The left 
plot shows the cases of (l/L) = 0.005 (o), 0.01 (T), and 0.02 (o) for a fixed resolution Air = 0.05. 
The right plot shows the results by the two different resolutions, Aa; = 0.05 (▼) and 0.025 (V), in 
the case (l/L) = 0.01. In the range x out < log(4L), the error value Sy does not depend on the 
resolution and decreases with the increase in x ou t- In the range x ut ^ l°g(4£), the error value 
8y depends on the resolution and becomes larger with the increase in x out in both resolutions. 
The error is interpreted as the boundary error and the nonsystematic error for x ou t l°g(4£) and 
x Q ut ^ log(4L), respectively. See text for details. 



coincide with those of the solution in [29|. For example, the figures in [29] indicate that 
Y > and Y takes the maximum value where the brane and the horizon cross each other, 
and the same feature holds in our numerical data. However, as (l/L) is increased for a 
fixed Xout) a strange behavior becomes relevant. Figure [3| shows our numerical solution 
of Y(x,x) for (l/L) = 0.06, x ou t = 5.0, and Ax = 0.05. The value of Y does not decay 
properly at the distant region. Moreover, there is an unnatural jump in the values of Y at 
the grids 3=1 and 2 (i.e. x = anci 2A%). This is in contrast to the case of a zero 
tension brane in Sec. |], where the error in Y is suppressed in the neighborhood of the axis. 
As the characteristic error in Y, we define the following quantity: 

5 Y := max \Y(x, 2 A X )- Y ext (x)\. (5.1) 

Here, l^ x t(x) means the extrapolated value of Y at x = 2A% using the values of y(x,0) 
and Y(x,Ax)- Let us observe how the value of 5y depends on x ou t and (l/L). 

Figure [| shows the value of 5y as a function of x out . In the left plot, the cases 
(l/L) = 0.005, 0.01, and 0.02 are shown for a fixed Ax — 0.05. As x ou t is increased, Sy 
decreases at first, but later increases. The change from decrease to increase happens at 
Xout — log(4L) (i.e. pout — 4£). This behavior is opposite to our expectation, since the 
error has to continue to decrease as x out is increased until it hits one of the relaxation and 
grid errors which are both < 10 -5 . Therefore, although the error 5y for x ou t ^ log(4L) 
can be interpreted as the boundary error, we have detected the nonsystematic error in the 
range x ou t ^ log(4L). The right plot shows the comparison between the two resolutions, 
Ax = 0.05 and 0.025, for a fixed (l/L) = 0.01. For x out < log(4L), the two error values 
almost coincide and they are determined only by x ou t- Therefore, we can confirm that 
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Figure 5: Dependence of the error Sy in Y on the values of (1/L). The left plot shows the cases 
of x out — 4.0 (x), 5.0 (A), and 6.0 (•) for a fixed resolution Ax = 0.05. The right plot shows the 
results by the two different resolutions, Ax = 0.05 (A) and 0.025 (A), in the case x ou t = 5.0. The 
boundary error 5y in the range (1/L) < 4exp(— x out ) scarcely depends on the resolution and grows 
linearly with respect to (1/L). The nonsystematic error in the range (1/L) > 4exp(— x ou t) depends 
on the resolution and grows nonlinearly with respect to (1/L) in both resolutions. 



the error in this region is the boundary error. On the other hand, for x ou t ^ log(4L), the 
error 5y depends on the resolution. However, they are not the grid errors, since the grid 
errors are < 10 -5 and ~ 10~ 7 for Ax = 0.05 and 0.025, respectively, as proved in Sec. |4|. 
The relaxation error is also ~ 10~ 6 . The important feature is that in both resolutions, the 
errors become larger as x ou t is increased. Therefore, we again confirm that the error in the 
range x ou t ^ log(4L) is the nonsystematic error. 

Figure |B| shows the value of 5y as a function of (1/L). The left plot shows the cases 
of x on t = 4.0, 5.0, and 6.0 for a fixed Ax = 0.05. We see the boundary error in the range 
(1/L) < 4exp(— x ut), arid it grows linearly with respect to (1/L) for a fixed x ou t- This is 
natural because the value of Y itself grows linearly with (1/L). On the other hand, we see 
the nonsystematic error in the range (1/L) > 4exp(— x ou t). In contrast to the boundary 
error, the nonsystematic error grows nonlinearly with respect to (1/L) for a fixed x ou t; arid 
the growth rate is roughly given by ~ (1/L) 4 . The right plot shows the comparison between 
the two different resolutions, Ax = 0.05 and 0.025, for a fixed x ou t = 5.0. Again, we confirm 
that the boundary errors are almost same for both resolutions, while the nonsystematic 
error depends on the resolution. Although the nonsystematic error for Ax = 0.025 is 
smaller than that for Ax = 0.05 if compared with the same (1/L), the nonlinear growth 
~ (1/L) 4 is observed in both resolutions. 



5.2 The error in the surface gravity k 

As another example of the numerical error, we look at the error in the surface gravity k. 
As mentioned in Sec. 2.2 , there are the four boundary conditions on the horizon, and the 
three of them are chosen in our numerical calculation. Since we omitted the boundary 
condition that corresponds to the constancy of k on the horizon, we have to check if it 
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Figure 6: The same as Fig. [| but for the error S K in the surface gravity. The boundary error and 
the nonsystematic error can be seen for x ou t ^ log(4L) and x ou t ^ log(4Z), respectively, also in the 
error S K . 
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Figure 7: The same as Fig. [| but for the error 5 K in the surface gravity. The nonlinear growth of 
the nonsystematic error with respect to (1/L) is seen also for S K . 

is satisfied within the amount of the systematic errors after generating the solutions. In 
order to evaluate the error in k, we define 

5 K := max[«] — min[ft]. (5-2) 

Figure || shows the value of 5 K as a function of x ou t- The left plot shows the cases of 
(1/L) = 0.005, 0.01, and 0.02 for a fixed Ax = 0.05, and the right plot shows the comparison 
between the two different resolutions, Ax = 0.025 and 0.05, for a fixed (1/L) = 0.01. 
Although 5 K is smaller than 6y, they have common features. Namely, if x ou t is increased, 
the value of 5 K decreases at first but begins to increase at x ou t — log(4L). Therefore, also 
for S K , we find the boundary error and the nonsystematic error in the ranges x ou t ^ log(4L) 
and x ou t ^ log(4L), respectively. Although the nonsystematic errors for the resolutions 
Ax = 0.05 and 0.025 have a little different values, they both become larger as x ou t is 
increased. 
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Figure shows the value of S K as a function of (l/L). The left plot shows the cases 
of x ou t = 4.0, 5.0, and 6.0 for a fixed Ax = 0.05, and the right plot shows the comparison 
between the two different resolutions, Ax = 0.05 and 0.025, for a fixed x ou t = 5.0. The 
boundary error for 5 K is almost constant when (l/L) is changed and this is different from 
the behavior of 5y oc (l/L). This is because the surface gravity k itself is almost constant. 
On the other hand, for a fixed x ou t; the nonsystematic error for 8 K grows nonlinear ly with 
respect to (l/L). Although the nonsystematic error for Ax = 0.025 is smaller than that 
for Ax = 0.05, we observe the nonlinear growth with respect to (l/L) in both resolutions, 
and again the growth rate is roughly ~ (l/L) 4 . 

5.3 Comparison with the perturbative study 

In the next subsection, we discuss the fact that the nonsystematic error has a crucial 
meaning for the existence of the solution. But in the perturbative level, the (at least C 1 ) 



solution of a black hole on a brane was found for a small mass [30, 31]. In this subsection, 
let us examine whether our numerical solution is consistent with the perturbative study 
if the numerical errors are ignored. Since the perturbative study |p0| 31] was done in 



a different gauge condition, the direct comparison between the metric functions is not 
possible. However, we can compare the two results by calculating nondimensional quantities 
on the horizon. 

As the first example, let us consider the product of the surface gravity k and the proper 
radius of the horizon on the brane 



p h := p^exp [C(p h ,w/2)] 



(5.3) 



(i.e., the proper circumference divided by 2tt). Here, C is written as a function of the 
coordinates (p,x)- Figure || shows the relation between (l/L) and K/5 D for fixed values of 
Ax = 0.05 and x ou t = 5.0. The prediction by the perturbative study l3l|] is shown by 
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a solid line (see Appendix A for a derivation). The two results agree well in the regime 
(1/L) « 1. 

As another example, let us consider the quantity related to the distortion of the horizon. 
For this purpose, we define the effective radius of the horizon in the 5-dimensional sense as 

p cS := (A 5 /n 3 f 3 , (5.4) 

and calculate the ratio of the effective radius p e R to the proper radius on the brane pb- 
Here, ^3 denotes the 3-dimensional area of a unit sphere, £1% = 2tt 2 , and A$ is the area of 
the horizon 

A 5 :=20 2 / - e R+2C sm 2 X dx, (5-5) 







with := 4-7T. Figure^ shows the relation between (1/L) and p e s/Pb for fixed values of 
Ax = 0.05 and x out = 5.0. The value of p e s/pb decreases as (1/L) is increased, and thus 
the horizon becomes flattened. The prediction by the perturbative study is shown by 
the solid line (see Appendix A for a derivation), and the numerical data agrees well with 
the perturbative study in the regime (1/L) <C 1. 



As found above, our numerical data is consistent with the perturbative study |3C], |31] 
if the numerical errors are ignored. This result can be interpreted as the evidence for the 
correctness of our code as well as the correctness of the perturbative study. 

5.4 Interpretation 



Now, let us discuss the interpretation of our results. In Sees. 154] and p.2[ we have found 
the nonsystematic errors in our numerical data for both the error Sy in Y and the error 
5 K in the surface gravity. These nonsystematic errors become relevant for x out > log(4L), 
and grow nonlinearly with respect to (1/L) for a fixed x ou t- On the other hand, as shown 



in Sec. 5^3, our results are consistent with the perturbative study |30], 31] for a small black 
hole on a brane in the regime (1/L) <C 1 if the numerical errors are ignored. In order to 
interpret these strange results, we suggest and discuss the following three possibilities: (I) 
The numerical code is wrong; (II) The solution does not exist; and (III) The solution exists 
but it is numerically unstable. 

We understand that the possibility (I) is the case that most readers might expect. 
However, we have successfully reproduced the Schwarzschild solution for (1/L) = in 
Sec. |j. Moreover, our results for (1/L) <C 1 are consistent with the perturbative study 



|3C, 31 1 . This guarantees that at least the linear terms in p/i are correctly included in our 
code. Therefore, our task is to properly include the other nonlinear terms, and this can be 
done by careful checks. For this reason, we definitely exclude the possibility (I). 

We turn to the possibility (II) that no solution exists. Here, we consider the case that 
no regular solution exists but the solution to the finite difference equations exist. Since 
such a case can happen, it is generally warned that at least the appropriate convergence 
has to be checked in numerical calculations. Let us recall Figs. || |6|, and 0. In the right 
plots of those figures, we compared the errors of the results by the two different resolutions, 
Ax = 0.05 and 0.025, and found that the numerical data does not show the appropriate 
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convergence. If we assume that the solution to the finite difference equations is unique (see 
possibility (III) for the meaning of this assumption), this indicates that a regular solution 
does not exist at least for finite x ou t- If the error continued to decrease with the increase 
in x ou t, the solution would exist for x out = oo. However, the fact is that the error becomes 
larger with the increase in x ou t in the range x ou t ^ log(4L), raising doubts for the existence 
of the solution also for x ou t = oo- 

For small (1/L), the solution exists at the perturbative level [31] and our numerical 
data is consistent with this study. But we also found the nonsystematic error that grows 
nonlinearly with respect to (1/L). Here, let us recall the fact that the regularity of the 
horizon is not perfectly guaranteed in the perturbative study (see also Appendix A). If 
the perturbation becomes singular at the higher-order perturbation, what happens in the 
numerical calculation? Such a singular effect should appear as the numerical error that 
behaves in a nonsystematic way. The nonlinear growth of the nonsystematic error with 
respect to (1/L) supports this view. We consider that this is the most natural interpretation 
of our result. If this is the case, the regular solution does not exist at least to this setup. 

It is interesting to ask at which order the perturbation becomes singular. This is a 
difficult question to answer, and we just make a speculation here. In Figs. |5| and [7|, we 
found that the nonsystematic error grows quite rapidly, ~ (1/L) 4 , as (1/L) is increased. 
If we assume that this power of the nonsystematic error reflects the order at which the 
perturbation becomes singular, the singular behavior might appear at the 4th-order per- 
turbation. 

Readers might wonder why we do not claim the existence of the solution at least for 
the range (1/L) <C 1 in spite of the consistency of our numerical data with the perturbative 
study and the small values of the nonsystematic error. In fact, we have reproduced the 
results of the perturbative study for (1/L) <C 1 with a typical numerical error < 0.1%. 
However, the merits in the numerical calculation is the ability to explore the nonlinear 
regime in (1/L), and our results obviously show that something is wrong at the nonlinear 
level. In such a situation, it is much more dangerous to claim the existence of the solution 
however small the detected error is. We consider that this is a common sense for numerical 
researchers who are serious in analyzing the numerical errors. 

Readers might also expect that the solution exists for small (1/L) and it vanishes at 
some critical point as (1/L) is increased, as suggested in Ref. [17]]. We do not consider 
that this is the case, because as shown in Figs. |] and [6|, the nonsystematic error becomes 
relevant for x ou t > log(4L) for very small (1/L) values such as (1/L) = 0.005. To support 
this view, let us look at the existing example by the author and his collaborator [36] 
where disappearance of a solution was actually observed. In that paper, we numerically 
studied the apparent horizon formation in high-energy particle collisions, and found that 
the solution of the apparent horizon exists for a small impact parameter b, but it disappears 
at the maximal impact parameter 6 max - As shown in Fig. 2 of Ref. fl36|] , there is a clear 
signal of 6 max such that dr m \ n /db diverges. Although the numerical error grows as b is 
increased, the error turns out to be the grid error whose growth is driven by the distortion 
of the apparent horizon. This is the typical phenomena that happens when a solution 
vanishes as the system parameter is changed. However, in the present setup, we do not 
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find any signal of disappearance of the solution. Besides, the numerical error that grows 
with (1/L) is not the grid error, but has turned out to be the nonsystematic error. These 
features are all in contrast to the case of |3(| . 

Finally, let us discuss the third possibility (III) that the solution exists but it is nu- 
merically unstable. Here, we refer to the case that a regular solution exists and that two 
or more solutions to the finite difference equations exist. Among these numerical solutions, 
only one solution approximates the regular solution, and our obtained solution is one of the 
other artificial solutions that do not represent the true solution. We consider that this case 
is unlikely, since our solution reproduces the result of the perturbative study for (1/1) < 1 
and hence it does not seem to be artificial. However, at the same time, we cannot rigorously 
show the uniqueness of the solution to the finite difference equations since they are highly 
complicated. As a result, we cannot rigorously claim the nonexistence of the solution, i.e., 
that the possibility (II) is the case. For this reason, our results also have to be carefully 
interpreted as well as the other numerical or perturbative studies. In order to show the 
nonexistence of solutions rigorously, a proof by the method of the mathematical relativity 
will be required. 



6. Summary and discussion 

In this paper, we studied a black hole on a brane in the RS II scenario using the nu- 
merical code having almost 4th-order accuracy. Although we successfully reproduced the 
Schwarzschild solution in the case of a zero tension brane (Sec. Eh, the nonsystematic error 



was detected for the cases of a nonzero tension brane (Sees. 5.1 and |5.2[ ). This nonsys- 
tematic error becomes relevant for x out ^ log(4L) even for a very small (1/L), and it 
grows nonlinearly with respect to (1/L) for a fixed x OU f On the other hand, our numerical 
data agrees very well with the perturbative study J30|, |3l| if the numerical error is ignored 



(Sec. |5.3|) . In Sec. p.4| , we discussed the interpretation of our result. The most natural 
interpretation is that the perturbation becomes singular at the higher order in (1/L), and 
this effect appears as the nonsystematic error in the numerical data [the possibility (II)]. 
If this is the case, no solution of a black hole on a brane exists. We also pointed out the 
remaining possibility that the solution exists [the possibility (III)], although it is unlikely. 
To summarize, our numerical results raise strong doubts for the existence of a static black 
hole on a brane. To be precise, our claim is: A solution sequence of a static black hole 
on an asymptotically flat brane that is reduced to the Schwarzschild black hole in the zero 
tension limit is unlikely to exist. It is also worth pointing out that a similar claim has been 
made in the context of the AdS/CFT correspondence 14]. 

In this paper, we focused our attention to the numerical calculations. So far, we have 
not succeeded in proving the nonexistence of solutions by an analytic method. But, here, we 
would like to give an analytic discussion that supports it. Let us ask in this way: Suppose 
the functions T, R, and C satisfy Eqs fl2.5|) - (l2.7D and fl2.9|) and the boundary conditions 
at infinity ( 2.11[) , at the symmetry axis (pTTH) (|2TT4|) , on the horizon ([D]|-(|l|), and 



T x /T = R x = (p/l)(e R - 1) on the brane. Then, can the function C satisfy the boundary 
condition C iX = (p/£)(e R — 1) on the brane? To answer this question, we calculate the 
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difference of the equations (|2.7| ) and (|2.S|), The result is the equation of the form 

~P 2 C,xx + C iPP + F[T XX , T x , R x , C tX , T p , R p , C iP , T, R, C,p, x] = 0. (6.1) 

Here, there is no singularity in this equation at x = 0. This is a wave equation where 
X is the "time" coordinate. The "initial condition" is given by C = R and C tX = at 
X = 0. For given T and R, we can "evolve" C using this equation from x = to tt/2. 
Then, the value of C x at x = 11 / 2 is determined as a result of the "temporal evolution", 
and therefore, there is no guarantee that C satisfies the junction condition on the brane 
C :X = (p/£)(e R — 1). Although this condition can be satisfied in a special situation (e.g. the 
spherically symmetric case), such a chance cannot be expected in general. We consider that 
this discussion is consistent with our numerical results. In Fig. ||, we found an unnatural 
jump in the value of Y in the neighborhood of the symmetry axis. This would be because 
the boundary conditions on the brane and at the symmetry axis are incompatible. 

Suppose solutions of a static black hole on a brane do not exist in the RS scenarios 
as discussed in this paper. Then, what happens after the gravitational collapse? As 
an example, let us consider a head-on collision of high-energy particles. Because of the 
nonexistence of the solution, one might expect that the black hole does not form in the 
collision. However, we do not consider that this is the case, since there exist studies on 



the apparent horizon formation in the RS II scenario [37, 38], and even the formation 
of a very large apparent horizon is reported |38|]. For this reason, the black hole should 
form. Although the produced black hole relaxes to a quasi-equilibrium state by emitting 
the gravitational radiation, it does not completely become static and cannot remain on 
the brane. In this sense, the black hole is unstable. We suggest two possibilities of the 
consequences of the instability: the event horizon pinch and the brane pinch. 

In the first possibility of the event horizon pinch, we refer to the case where the brane 
scarcely moves and the cross section of the event horizon and the brane shrinks. In the 
classical general relativity, the pinch of an event horizon is forbidden (e.g. [|35[|). However, 
if the radius of the cross section becomes order of the Planck length, the event horizon 
pinch could happen by quantum gravity effects and the black hole may leave the brane. 
We consider that this scenario may happen when the horizon radius is sufficiently larger 
than the bulk curvature scale, since the brane is expected to be rigid in this case. 

In the second possibility of the brane pinch, we refer to the case where the event 
horizon scarcely moves while the brane moves. The motion of the brane continues until 
it pinches and generates a baby brane. We consider that this scenario may happen when 
the horizon radius is smaller than the bulk curvature scale, since in this case the brane 
can easily fluctuate. Basically, this scenario was already suggested and studied in detail 



in HO, 41, HJ| by numerical simulations for a zero tension brane and a more detailed 
model of a brane (but without the self-gravity effects). In those studies, the authors gave 
the brane an initial velocity in the bulk direction, and observed the brane pinch. The 
new insight from the results in the present paper is that such a brane pinch might happen 
because of the instability without the initial velocity in the bulk direction. We are also able 
to estimate the time scale of this phenomena. In our numerical data, the nonsystematic 
error has the order of ~ (1/L) 4 . Since this error amount is expected to cause the temporal 
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evolution, the typical time scale would be ~ £ 4 / 1 p\. Therefore, the brane pinch may happen 
relatively slowly unless the velocity in the bulk direction is given. In order to clarify if this 
is the case, a simulation taking account of the self-gravity effects of the brane is necessary. 
It also has to be checked if the brane pinch can happen when the Z^-symmetry is imposed. 

The above discussion has the following implication for the black hole production at the 
LHC. If the extra dimension is warped, the produced black hole may escape into the bulk. 
If this is the case, it is difficult to observe the signals of the Hawking radiation especially 
from black holes with large mass, since the time scale of the instability may be shorter 
than the time scale of the evaporation. Therefore, the instability of black holes on a brane 
gives us many interesting and important issues to be explored. 
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A. A brief review of the perturbative study 



Karasik et al. [30, 31 1 solved the static solution of a black hole on a brane in the RS 
scenarios in the case of a small mass by a perturbative study. Here, we summarize the 
part of their results that is directly related to the discussions in this paper. In this section, 
we follow their notation. Since they were mainly motivated by the TeV gravity scenarios, 
they studied in the framework of the RS I scenario. However, their results can be applied 
also to the RS II scenario by changing the sign of I. 

They considered the situation where the radius of the black hole p := ^SG^M/Sir is 
much smaller than the bulk curvature scale £, and adopted the so-called matching method. 
In this method, the spacetime is a perturbed Schwarzschild black hole in the neighborhood 
of the horizon and is a perturbed AdS spacetime at the distant region, where e := pjt is 
the small parameter. In the matching region p, <C r -C I, the two approximations both 
hold and the two solutions can be matched to each other. 

In the neighborhood of the horizon, the metric is assumed to be 

2 



A 
B' 



-Bdt z + ^dp 2 + 2Vdpdip + p 2 Udip 2 + p 2 sin 2 ipd^l 



(A.l) 



where z := i — p cos ip and A, B, V, U are the functions of p and ip. The background 
solution is the Schwarzschild solution Aq = Uq = 1, Bq = 1 — 1/ p 2 , and Vq = 0. The first 
order quantities for these functions are given as eA\, — ef?i, eU\ and eV±, respectively, and 
it turns out that they are expressed in terms of F, H, F „, H tP , F^ and -ff »/>, where F is a 
gauge function and H is a wave function. The wave function H satisfies the equation 

(p 2 - 1) ( H pp + 2 C °*y + 1 l Hi, = 0. (A.2) 



p J ' sin V> cos ip 
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It was found that when the perturbed AdS metric at the distant region includes only 
integer powers of M, the matched solution in the neighborhood of the horizon is singular 
at p = 1. However, if we allow the terms with half integer powers of M at the distant 
region, the metric becomes at least C 1 and the zeroth and first law of the black hole 
thermodynamics can be satisfied. The solution of H is 

h = 5 - Aq w + q {p) (JL g m^) + _ * (3) W 

07T ' \ 6TT Sill 1p O 



+ 3 r^+^^r 



5 ( 2 )(V) +^a(A)Q (A _ 1)/2 (p)G(A,V). (A.3) 

n=lJ A=3 



Here, P„ := p v / ^lP I/ 1 (2p 2 - 1), and & := py^lQi (2p 2 - 1), and <? (2) , S (3) and 
G(\,ip) are some elementary functions of ip (see Eq. (36) in (3^] for details). The values of 
do, d,2, and a(A) are fixed as 5q = d?, = and a(A) = — (16/37r)(— 1) a /(A 2 — 4) by imposing 
boundary conditions and the zeroth and first laws of the black hole thermodynamics. The 
expansion of the function Q u has the terms proportional to (p 2 — l) n log(/9 2 — 1) for integer 
n > 1. Hence, H^ p has the form 

# P ^ ^7 [log(p 2 - 1) + 2 1E + 2^(1/2)] + fnty), (A.4) 
sin ip 

in the neighborhood of the horizon, where ho(ip) and h\{ip) are given by Eqs. (44) and 



(45) in Ref. [31]. Here, a miraculous cancellation happens and ho(ip) becomes zero. This 
means that the term proportional to (p 2 — 1) log(/) 2 — 1) is not included in H. Then, the 



equation ( A. 2 ) implies that H pp is also finite. As a result, the function H includes terms 
proportional to (p 2 — l) n log(/) 2 — 1) only for n > 3. 

When the horizon is regular, the location of the horizon pn and the surface gravity n 
are given by 

p H = 1 + e [-6H(1, tt/2) + H p (l, tt/2)] ; (A.5) 
K=l + 6eH(l,ir/2), (A.6) 



from Eqs. (23c) and (23d) in [31], and the 3-dimensional area of the horizon is given by 

^ = l-|cff(l,7r/2) (A.7) 



from Eqs. (27) and (30) in [31], where Ajj := 2ir 2 p 3 . From these formulas, we immediately 
find that the product of the surface gravity and the proper radius of the horizon on the 
brane is npn = 1 + eh\{ir/2). On the other hand, the ratio of the effective horizon radius 
to the radius on the brane becomes (A H /A^) 1 / 3 / p H = 1 + e [(9/2)^(1,^/2) - /ii(vr/2)]. 
Here, H(l,ir/2) = 8/9ir from Eq. (38) in |3l| and hi(ir/2) is numerically evaluated to be 
~ 1.087336. In order to apply to the RS II scenario, we change the sign of I and find 

Kp H * 1 - 2.17467 x V2G 5 M/3* . (A g) 



- 22 - 



(Ohl/3 



1 - 0.371801 x 



(A.9) 



PH i 
These formulas are used in Figs. || and ||, respectively. 

We comment on the remaining issues in this study. In order to guarantee the regularity 

on the horizon, the Kretchmann scalar K := R pupa R^ upcT = + eif W + e 2 K^ H 

has to be finite on the horizon. If the function H includes a term proportional to (p 2 — 
l) 3 log(/> 2 — 1), the third derivative H jPpp logarithmically diverges at p = 1. Although the 
authors of Ref. have confirmed the absence of the term H ppp in K^ 1 ' , there remains 
a possibility that diverges. In order to check this, two things are required. First, 
one needs to confirm if H )Ppp really possesses such a logarithmic term. But this is difficult 

has the infinite summation of 



because the coefficient of (p — 1) \og(p — 1) in Eq. 
terms whose absolute values become unboundedly large as A is increased. It turns out that 
this infinite summation does not converge even in the sense of distribution. Therefore, 



there is an ambiguity in the interpretation of Eq. (A.5). In the particle physics, several 
procedures to regularize diverging infinite summations are known. Although it may be 
possible to regularize the infinite summation in H in a similar way, this is not an easy task, 
and this situation makes it difficult to confirm the (non)existence of the term proportional 
to (p 2 - l) 3 log(/) 2 - 1) even numerically. Next, the presence of a logarithmic term in 
H ppp does not immediately imply the divergence of the Kretchmann scalar K, since the 
terms from the 2nd-order perturbations may cancel this term. However, the study of the 
2nd-order perturbation is also a difficult task. 

Finally, even if the regularity of the perturbation is established to some order, there 
remains the possibility that the metric becomes singular at the next order. If our in- 
terpretation in Sec. |5,4| is correct, the singular behavior might appear at the 4th-order 
perturbation. Therefore, it may be very difficult to check the (non) existence of solutions 
of a black hole on a brane by the perturbative method. 
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